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III.  Summary  of  Results 

In  previous  documents,  especially  the  Semi-Annual  Technical 
Reports  No.  2 and  3 submitted  in  June,  1974,  investigators  on  the 
ARPA  Climate  Dynamics  Project  at  Brown  have  reported  in  detail 
their  progress  towards  our  stated  objectives.  In  this  Final  Re- 
port, it  seems  approrriate  to  review  our  major  accomplishments, 
and  to  stress  the  general  impact  of  our  wcrk,  beyond  the  indiv- 
idual papers  cited  in  earlier  reports  and  listed  below.  Six 
major  accomplishments  can  be  cited: 

1.  Completion  of  the  first  ARPA-CLIMAP  18,000  YBP  numerical 
climate-simulation  experiment.  This  experiment,  the  concept  and 
design  of  which  originated  in  discussions  among  ARPA  investiga- 
tors at  Brown  and  at  RAND,  was  carried  out  in  a unique  enterprise 
which  involved  more  than  100  people;  a dozen  or  more  institutions; 
and  fiscal  support  both  by  ARPA  and  by  the  IDOE/CLIMAP  project. 
ARPA's  contribution  at  Brown  was  the  funding  of  part  of  the  work 
of  the  three  principal  investigators,  whose  efforts  were  central 
to  the  experiment  Its  contribution  through  RAND  involved  the 
support  of  the  RAND  General  Circulation  Model. 

This  experiment  marks  a turning  point  in  the  study  of  the 
dynamics  of  long-term  climatic  change.  For  the  first  time  a 
substantial  data  base  was  gathered  (mainly  under  the  CLIMAP  pro- 
ject) to  fix  the  sea-surface  temperature  pattern  during  an  ice 
age.  This  input,  together  with  ice-margin  and  sea-level  data. 
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form  the  boundary  conditions  for  the  running  of  an  atmospheric 
model.  The  model  then  calculates  many  elements  of  the  atmos- 
pheric circulation  which  maintained  the  ice-age  climate.  The 
experimental  aspect  of  the  project  lies  in  the  fact  that  ARPA 
scientists  at  Brown  gathered  independent  data  on  the  actual 
18,000  BP  response  of  the  atmosphere  at  point  locations  over 
the  continent.  By  comparing  the  model  calculations  with  inde- 
pendent geological  observations,  we  can  in  principle  confirm 
or  deny  the  accuracy  of  the  atmospheric  model  — and  thus  get 
at  two  fundamental  problems:  verification  of  the  explanatory 

power  of  numerical  models  of  the  atmosphere;  and  discovery  of 
atmospheric  mechanisms  which  maintain  our  planet's  climate  in 
other  modes  than  are  known  from  the  past  century  of  instrumental 
observations. 

Thus,  in  our  opinion,  the  results  of  the  Brown  ARPA  Cli- 
mate Dynamics  project  have  a significance  beyond  that  which  might 
appear  in  a simple  listing  of  our  publications.  At  a conference 
held  at  Brown  on  March  26,  1975,  25  scientists  gathered  to  ob- 
serve the  confrontation  of  geological  data  with  the  calculations 
of  the  RAND  model.  The  main  results  were  three:  First,  the 

surface  temperature  calculations  of  the  model  appear  to  be  con- 
firmed by  geological  observations  (See  Appendix  A) . Second, 
the  average  air-temperature  change  over  the  lands  during  the 
18,000  BP  ice-age  N.  Hemisphere  summer  was  about  -6.5°C.  The 
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importance  of  this  figure  lies  in  its  small  magnitude.  Consider- 
ing that  the  Northern  Hemisphere  during  the  last  century  has 
oscillated  through  a range  of  l.G°C  (annual  mean  values),  the 
surprise  is  that  an  ice-age  anomaly  should  be  so  small.  The 
general  conclusion  to  be  made  is  that  the  Earth's  climate  should 
be  considered  in  rather  sensitive  balance.  Third,  the  model 


precipitation  anomaly,  at  least  in  its  geographic  pattern,  is 
not  confirmed  by  geological  observation. 

2.  The  principal  investigator  of  the  Drown  ARPA  Project 
participated  during  June  1974  in  a two-week  session,  held  in  the 
USSR,  of  Working  Group  VIII  of  the  US-USSR  Bilateral  Agreement 
on  the  Environment.  As  a result  of  these  negotiations,  a bi- 
lateral program  of  paleocl imatic  data  exchange  was  initiated 
(See  Appendix  B; . The  importance  of  this  lies  in  the  need  to 
obtain  (for  the  running  and  verification  of  numerical  simula- 
tions of  past  climates)  better  data  on  the  ice-margins  and  pol- 
len records  from  the  vast  interior  lands  of  Russia.  Sometime 
in  1975,  or  possibly  early  1976,  a working  conference  of  paleo- 
climatologists  is  to  be  held  in  Moscow.  Hopefully,  that  con- 
ference (which  is  now  being  finalized  by  negotiations  of  NOAA 
representatives)  will  actually  materialize.  In  the  meantime, 
contacts  with  M.  G.  Grosswald  at  the  Institute  of  Geography  in 
Moscow  have  already  borne  fruit  in  the  form  of  more  accurate 
information  on  the  position  of  ice-margins  in  Siberia  during 
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i ; 

the  10,000  BP  ice  age. 

3.  The  i irst  units  of  a computer  data  bank,  in  which  in- 
formation on  moder,.  and  fossil  pollen  are  stored,  have  been  com- 
pleted. This  bank  has  already  been  drawn  upon  for  two  kinds  of 
analysis  (See  Appendix  C) . First,  the  writing  of  transfer  func- 
tions which  make  possible  the  extraction  from  pollen  data  quanti- 
tative estimates  of  past  climates.  Second,  the  preparation  for 
publication  of  several  synoptic  studies  of  climate-related  changes 
in  vegetation.  Up  to  now,  most  pollen-based  studies  of  past  cli- 
mates have  focussed  on  down-core  examination  of  single  sites. 

The  long-range  significance  of  our  program  is  its  emphasis  on 
the  discovery  and  documentation  of  synoptic  patterns  of  clima- 
tic change . 

4.  Our  technique  of  studying  simultaneously  pollen  grains 
and  planktonic  skeletons  in  the  same  samples  of  deep-sea  cores 
has  opened  the  way  for  a better  look  at  the  dynamics  of  clima- 
tic change.  In  a test  study  off  Cape  Hatteras,  deep-sea  cores 
were  used  in  this  way  to  obtain  evidence  on  the  timing  of  vege- 
tation changes  on  land  compared  to  the  timing  of  changes  in  the 
position  of  the  Gulf  Stream  (See  Appendix  D) . Evidence  is  that 
the  oceanic  changes  here  precede  changes  in  the  adjacent  contin- 
ent. Thus  encouraged,  ARPA  investigators  at  Brown,  in  coopera- 
tion with  palynologists  at  NYU,  began  a similar  study  in  cores 
eff  the  coast  of  California,  Oregon,  and  Washington.  Here  con- 
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tmuous  vegetation  records  can  be  obtained  going  back  about 
100,000  years.  Again,  we  are  confident  that  this  work  will 
prove  to  be  a springboard  for  much  future  research.  On  April 
25,  1975,  a conference  of  22  investigators  was  held  at  Brown  to 
consider  the  opportunities,  assess  results  to  date,  and  encour- 
age others  to  adopt  the  procedure. 

5.  Investigators  on  the  Brown  ARPA  Project  have  played  a 
leading  role  in  development  of  transfer-function  techniques  by 
which  geological  observations  on  pollen  or  marine  microfossils 
in  ancient  sediments  can  be  translated  into  quantitative  esti- 
mates of  climatic  variables.  Early  papers  by  Imbrie  and  Kipp 
(1971)  and  Webb  and  Bryson  (1972) , which  antedated  our  Project, 
initiated  a wide  and  growing  effort  by  many  investigators  to 
employ  these  techniques  to  document  the  climatic  history  of  the 
Quaternary.  With  ARPA  support,  we  have  carried  on  a systematic 
analysis  of  the  strengths  and  weaknesses  of  transfer  function 
methodology,  examined  basic  assumptions,  and  carried  out  alge- 
braic experiments  to  test  the  attributes  of  alternate  techni- 
ques. One  summary  paper  (See  Appendix  E)  is  in  press.  Another 
is  in  preparation  by  the  same  authors. 

6.  Climatic  history  without  an  accurate  absolute  chronology 
is  of  little  value.  Absolute  dates  based  on  Thorium-growth  methods 
from  fossil  reefs  on  oceanic  islands  have  proved  to  be  the  key- 
stone block  of  our  knowledge  of  the  chronology  of  the  climatic 
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history  o.l  the  past  250,000  years.  Again,  Brown  investigators 
played  a major  role  in  this  area  before  the  ARPA  Project  here 
was  funded.  During  the  life  cf  our  Project  here,  attention  was 
focussed,  first,  on  improving  the  accuracy  of  the  chronology; 
and,  second,  on  extracting  information  on  the  rate  of  change  of 
sea  level,  hence  on  the  rate  of  change  of  polar  ice  caps.  Major 
results  on  the  first  objectives  have  been  published  (notably  in 
Bloom  _et_  a 1 . , 1974)  and  been  reviewed  in  previous  ARPA  Reports. 

As  the  Project  terminated,  eight  new  Thorium  dates  led  to  the 
conclusion  that  the  60,000-year  high  stand  of  the  sea  is  not 
represented  in  the  Barbados  sequence.  And  the  thrust  of  our  work 
now  shifts  to  the  laborious  task  of  improving  the  accuracy  of 
the  Th-230  dating  method.  Major  results  of  the  second  objective, 
p-.eviously  reported,  are  that  the  125,000-year  sea-level  high 
stand  lasted  not  longer  than  5,000  years;  and  was  followed  by 
a sea  level  lowering  at  a rate  of  5-10  meters  per  thousand  years. 
The  dynamics  of  this  portion  of  the  sea-level  record  are  of  par- 
ticular interest,  because  the  climate  today  is  in  many  respects 
a duplication  of  planetary  conditions  obtaining  at  125,000  YBP. 
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PRELIMINARY  REPORT  OF  THE  CONFERENCE  HELD  AT  BROWN  UNIVERSITY 

MARCH  26,  1975 

i.  preliminary  conclusions 

A.  Abou t model  skill . 

1.  Model  surface-level  air-temperature  anomalies  over 
land  are  in  accord  with  generalized  geological  expectations;  and 
in  the  limited  number  of  sites  where  we  can  make  spot  checks, 
model  values  are  consistent  with  proxy  information. 

2.  Model  precipitation— anomaly  patterns  and  surface- 
wind  maps  are  not  consistent  with  geological  information  in  a 
number  of  areas.  This  conclusion  is  in  accord  with  the  relatively 

low  level  of  model  skill  in  simulating  these  qualities  for  today's 
world . 

B.  About  improving  verification  procedures  for  the  next 

experiment . 

3.  First  priority  should  be  given  to  increasing  the 
accuracy  and  distribution  of  data  from  pollen,  marine  mineralogy, 
loess,  and  Indian  Ocean  sediments. 

4.  Second  priority  should  be  given  to  opportunities  for 
improving  inferences  from  periglacial  features,  and  for  disen- 
tangling temperature  from  precipitation  information  in  snow-line 
data. 


5.  Model  readout  should  include  moisture-balance 


maps. 
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as  this  quantity  better  matches  the  pollen  verification  data. 

6.  Wind-roses  at  least  for  selected  sites  south  of  the 
ice  sheets  should  be  included  in  model  output. 

7.  The  model  used  to  calculate  wind  maps  at  the  surface 
should  be  augmented. 

8.  Greater  attention  should  be  paid  by  the  geologists 
to  the  degree  to  which  a site  represents  a 500-km-size  area, 
for  that  is  the  spatial  scale  of  model  response. 

C.  About  the  18K  v?orld . 

9.  From  the  proxy  data  assembled  for  verification  we 
have  learned  from  mineralogy  that  the  trade  wind  belts  shift 
equatorwards  about  5°  latitude;  and  from  palynology  described 
several  spatial  patterns  of  moisture  anomaly. 

10.  From  the  combined  operation  of  model-input,  model- 
output,  and  verification,  we  are  developing  a solid  basis  for 

least  the  broad  outline  > of  the  surface— temperature 
anomaly  of  the  sea  surface  and  the  surface-air  temperature  ano- 
maly during  the  last  glacial  maximum.  From  this  experiment, 
e.g.  an  August  sea  anomaly  of  2 . 3°C  correspo  nds  to  an  atmospheric 
anomaly  over  the  ice-free  land  of  about  6°C . 

II*  ^ summary  of  the  validation  phase  of  the  experiment  is  pro- 
vided by  Tables  1 and  2 on  the  following  two  pages. 


! 
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A*  Geographic  distribution  of  sites  and  site-clusters  for 
whibh  climatic  estimates  for  10K  v/cre  compiled  by  Webb  and  Krause. 


Region 

Temperature 

Anomaly 

Moisture 

Anomaly 

Totals 

Continent 

Latitude*1 

Annual  July 

N.  Amer. 

H 

2 

7 

10 

19 

L 

1 

1 

4 

6 

S . Amer . 

H 

— 

— 

1 

1 

L 

1 

- 

3 

4 

Europe 

H 

5 

4 

5 

14 

Asia 

H 

5 

5 

_ 

10 

L 

1 

mm 

1 

Australia 


H 

L 


Africa 


H 

L 


1 

3 


1 

5 


2 

.8 


Totals 


20b  18° 


29< 


67 


(a)  H =>34°  in  N.  Hemisphere,  >26°  in  S . Hemisphere 
L » between  34°N  and  26°S 

(b)  12  sites  from  pollen  data 

(c)  6 sites  from  pollen  data 

(d)  mainly  pollen  data 
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TaVlc*  2.  Generalized  geologic  estimates  of  maximum-glacial  continen- 
tal climatic-anomaly  contained  in  four  well-known  summary  publica- 
tions, compared  with  RAND  model  calculations. 


I 


Geological  Estimates 
Summary  I Quantity  Estimated 


Pul-Iication 


Ref . 
page 


Model 

Result  Calculations 

(T.  decrease  (N.  Hem.  summer, 
0<-:)  S.  Hem.  winter) 


Flint 

(.1971) 


Mean  annual  T in  mari- 
time climates  at  low  4 

altitudes  and  latitude? 

Maximum  mean  annual 
temperature  in  con-  4 

tinental  situations 


3°  - 6° 


4°  - 8°  (a) 


or  more  14°  (g) 


average  annual  glo- 
bal mean  (T) 

Central  Europe  (T) 


prob.  not 
over  7° 


6.4  (b) 


10°  - 14° 


8°  - 140  (c) 


West 

(1968) 


Precipitation  re- 
duction in  British 
Highlands 

Mean  annual  T 

Global  July  T 


15%  (global 
average) 


5°  - 8° 


> 4 . 5°  (h) 


6.4°  (b) 
6.4°  (b) 


Emiliani 

-and-Flint 

(1963) 


Central  Europe  (T) 
~N . Italy  (T) 
Carpathian  Poland 


8°  - 10° 


8°  - 14°  (c) 

Io°  (d) 

~9°~  (e) 


Emiliani  • 
0371) 


Continental  temp,  in 
low  latitudes  (6°5- 
27°N) 


5°  - 10° 


5°  - 7°  (f) 


(a)  Eastern  Coast  of  S.  America  north  of  22°S . Range  along  Atlantic 
African  Coast  is  3°  - 9°,  with  the  higher  figures  at  upwelling  sites 
tb)  Average  over  both  Hemispheres  of  ice-free  land. 

Xc)  Over  Europe*  latitudes  46°  - 50°. 

(d)  Mean  of  two  points  representing  Italy. 

(e)  If  mtn  point  to  east  is  averaged  = 11.5° 

(f)  Range  of  averages  of  corresponding  latitude  circles. 

(g)  Maximum  N.  Hera,  value  away  from  high  mtns. 

• (h)  A minimum  estimate  calculated  from  an  average  snowline-lowering 

of  760  m;  lapse  rate  of  0.6°C  per  100  m;  and  the  assumption  that 
precipitation  did  not  change. 
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REPORT  OF  THE  FIRST  MEETING  OF  WORKING  GROUP  VIII 
(INFLUENCE  OF  ENVIRONMENT  CHANGE  ON  CLIMATE)  - 
FOR  DAFLEHENTAT  1 ON’  OF  THE  AGREEMENT  BET’.VEEN 
THE  UNITED  STATES  OF  AMERICA  AND  THE  ULl6N  Or  S0/1ET 
SOCIALIST  REPUBLICS  ON  COOPERATION  JN  TH:a  PROT]S'i].0N 

OF  THE  ENVIRONMENT 

i T>3'T'T>nniir,.n''i  otj 


1,1  The  first  mooting  of  Working  Group  VIII  was  held  in  Le 
ningrad  on  Jure  10-21,  1974,  an  agreed  to  at  the  November  16, 

1975  meeting  of  the  Joint  Committee.  0 c participants  an  this 
meeting  from  both  the  U.S.  and  the  USSR  are  listed  in  Appendix  A. 
Plenary  sessions  v/erc  held  on  June  10-12  and  June  19-21.  Details 
of  the  agenda,  including  subjects  and  authors  of  formal  presen- 
tations, are  given  in  Appendix  B. 

During  the  period  June  13-18, ' nenbers  of  the  American  dele- 
gation visited  Soviet  institutions  and  observatories  in  Moscow, 
Novosibirsk.  Kiev.  Crimea  and  the  northern  Caucasus.  There  were 


also  visits  to  a number  of  institutions  in  Leningrad  during  the 
first  period  of  the  meeting.  The  institutions  visited  by  the 
American  delegation  are  listed  in  Appendix  B. 


1.2.  Working  Group  VIII,  by  decision  of  the  Joint  Committee, 
io  divided  into  three  projects.  To  facilitate  discussion,  the 
projects  have  been  organized  into  sin:  subgroups  as  follows. 
Project  1.  Joint  studies  of  the  effects  of  changes  in  the  heat 
balance  of  the  atmosphere  on  climate. 

a.  Subgroup  1 - The  influence  of  changes  of  solar 

activity  on  the  climate. 

b.  Subgroup  2 - Modelling  of  the  climate. 

c.  Subgroup  3 - Assessment  of  past  climates  on  the 

basis  of  the  analysis  of  natural  objects  and 

data. 

Project  2.  Joint  studies  of  the  effects  of  pollution  of 

the  atmosphere  on  climate. 

a.  Subgroup  4 - Monitoring  atmospheric  constituents 

and  assessment  of  their  effects  on  climate. 

b.  Subgroup  5 - The  effect  of  pollution  of  the  upper 

levels  of  the  atmosphere  on  climate. 


I 
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Project  2 (aloo  Subgroup  6).  Joint  studies  o'?  the  meteorology  and 
aii^-cea  interaction  of  polar  regions  in  both  hemispheres 
. ac  they  affect  the  climate  of  the  planet. 

There  were  discussions  concerning  the  need  to  recommend  the 
reorganisation  of  this  Working  Group  to  the  Joint  Commission.  Im- 
portant considerations  for  the  organization  of  the  Working  Group 
are  achieving  sufficient  disciplinary  homogeneity  within  each  pro- 
ject to  pex-nit  successful  collaboration  end  the  recognition  of 
the  relationship  of  this  Working  Group  to  other  groups  established 
under  this  and  other  bilateral  and  international  agreements. 

For  the  latter  reason,  the  work  of  the  Subgroup  on  polar  ctu-  - 
dree  is  specifically  concerned  with  problems  of  climate  and  not 
with  the  broad  spectrum  of  all  polar  investigations.  On  the  other 
hand  there  are  aspects  of  air-sea  investigations  without  wnich 
full  understanding  of  the  environmental  influences  on  climate 
cannot  bo  adequately  treated. 

At  the  same  time,  the  studies  of  solar  activity  on  climatic 
variations  use  different  types  of  data  and  analysis  and  are 
broader  in  scope  than  the  work  of  the  other  subgroups  of  Project  1. 

For  these  reasons  a recommendation  is  made  to  the  Joint  Commit- 
tee for  reorganizing  Working  Group  VIII.  It  is  realized  that  this 
recommendation  impacts  the  interface  between  the  efforts  under 
this  bilateral  agreement  and  that  on  the  Study  of  the  World 
Oceans.  It  is  therefore  anticipated  that  action  on  this  recommen- 
dation would  require  communication  between  the  Joint  Committees 
under  the  two  Agreements, 

The  following  revised  structure  of  Working  Group  VTII  is 
recommended: 

Section  1.  Joint  studies. of  the  effects  of  chanses  of  the  heat 
balance  of  the  atmosphere  on  climate, 
a.  Subgroup  1 - Modelling  of  climate. 

. b*  Subgroup  2 Assessment  of  past  changes  of  the  cli- 
mate on  the  basis  of  analysis  of  natural  objects 
and  data. 

c.  Subgroup  J>  - Interactions  of  the  atmosphere  with 
polar  regions  and  the  oceans  as  they  effect 

climate. 
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c)  The  properties  of  the  climatic  models  and  the  li- 
mits of  their  applicability  may  differ  depending 
on  the  methods  used  for  parameterization  of  phycica.' 
processes  in  the  atmosphere  and  the  number  of  the 
. feedbacks  incorporated.  It  is  necessary,  therefore, 
to  formulate  the  models  compatibly  to  the  require- 
ments and  to  evaluate  the  limits  of  model  appli- 
cability with  respect  to  concrete  problems, 

f)  Consideration  should  be  given  to  the  possible  de- 
sign, performance,  analysis  and  interpretation 
of  numerical  experiments  requested  by  other  sub- 
groups in  order  to  assess  the  role  of  specific 
external  (nan-made)  factors  effecting  the  climate. 

2.2.2.  Data  Exchange 

Bearing  in  mind  that  ail  the  current  theories  of  cli- 
mate arc  approximate,  it  is  necessary  to  check  these  theories 
against  empirical  material  before  applying  the  models  to  cal-  n 
culations  of  future  climatic  condi cions.  For  this  purpose, 
short-term  seasonal  or  intcrannual  observations  of  meteorolo- 
gical elements  are  needed  as  well  as  the  data  on  contemporary 
climate  variations  and  on  climate  change  in  the  geological 
past.  The  specification  of  new  observational  requirements  arc 
also  needed. 


2,2,3  Actions  Proposed 

A joint  r.vnposiun  on  climatic  noddling  may  be  convened 
in  Uzbekistan,  USSR  in  1976. 

The  United  States  delegation  would  like  to  issue  invi- 
tations to  Soviet  scientists  to  visit  the  following  IJ.S.  insti- 
tutions for  periods  of  6 to  12  months.  The  invitation  is  of- 
fered to  the  indicated  scientist  or  a colleague  to  be  agreed 
on  by  both  countries, 

a)  To  come  to  the  Geophysical  Fluid  fgnau.ics  Labora- 
tory, National  Oceanic  and  Atmospheric  Administra- 
tion, Environmental  Research  Laboratories  in  F„ im- 
ceton,  N.J.  (one  or  two  of  the  visits  to  take 
place  in  1975) . 


t 


L 


i 
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1'i-0j . ossor  L.I.Budyko  - Climate  sensitivity  and  stability- 
Pcoiossor  K.  la. Kondratiev  - Radiation  effects  of  aerosols 
Academician  G.I.Lnrchuk  - Numerical  integration  methods 
Professor  S.S. Zilitinkevich  - Coupled  ocean-atmosphere  models 

b)  To  come  to  the  National  Center  for  Atmospheric 
Research  Boulder,  Colorado. 

Dr#  A.F.Ti-c  ~hnikov  — Effects  of  sea— ice 
/ cademiciar  G.I.Larchuk  — adjoint  methods  applied  to 

climatic  variations 

c)  To  come  to  the  University  of  Washington 
(Seattle,  Washington) 

* ■ • 

Dr#  A,F .Trechnikov  — Modelling  of  ice  dynamics 

d)  To  come  to  Stanford  University,  California 
Academician  Mustol  - Problems  of  solar  activity  and  climate 


Subfir°«TLJ^_Asncfisrient  of  past  climates  on  the  b sis 
•2£ — ihc  analysis  of  natural  objects  and  d a t a 


2.J.I.  Scientific  Cooperation 

A.  Century  program  (0  - 100  years  b6fore  present) 

This  v/ill  focus  on  climatic  changes  in  the  period  of  mo- 
dem instrumental  data,  including  aerological  data  for  the  past 
20  or  50  years,  and  will  attempt  to  document  the  most  recent 
trends  of  climate  together  with  environmental  changes  which  may 
•have  causal  connections  with  such  trends. 

a)  Exchange  of  data,  as  raw  (unsmoothed)  time  series: 

(i)  Indices  of  atmospheric  circulation  intensity 
and  patterns  in  the  Northern  Hemisphere. 

- Update  of  Dzerdzeevskii  circulation-ty£e 
epoch  data  by  the  Institute  of  Geography, 
Moscow;  (Ya._L._ Rauner,  suggested  coopera- 
ting scientist) . 

- Updated  series  of  coefficients  of  empiri- 
cal orthogonal  functions  of  monthly  mean 


1 


r 

\ 
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sea-level  pressure  (by  NOAA  Environmental 
Data  Servico,  Washington,  D.C.;  Mitchell, 
Jr.,  suggested,  cooperating  scientist). 

(ii)  Regional  and  hemispheric-scale  average  values 
of  monthly  average  climatic  parameters  in  Nor- 
them  Hemisphere. 

- Climatological  surface  air  tempo  rat 'ire  in. 
the  form  of  digitized  grid-point  values  on 

' I IBM-compatible  magnetic  tape  for  the  period 
1801  to  at  least  I960  together  with  related 
publications  (by  L'ain  Geophysical  Observato- 
ry, Leningrad;  Ye. P. Borisenkov,  suggested 
cooperating  scientist) 

- Updated  series  (since  195b)  of  hemispheric  . 
circulation  statistics,  including  zonal  and 
meridional  wind  indices,  free-air  tempera- 
ture, geopotential  height,  specific  humidi- 
ty, kinetic  energy,  fluxes  of  heat  and  mo- 
mentum (by  KOAA/CFDL,  Prince  ton,  U.  J. ; 

A. Oort,  suggested  cooperating  scientist). 

- Complete  digitized  "World  Weather  Records" 
temperature,  precipitation  and  pressure  da- 
ta, by  stations,  on  magnetic  tape  for  the 
period  from  beginning  of  record  at  each 
station  to  at  least  I960  together  with  re- 
lated publications  (by  NOAA  Environmental 
Data  Service,  Washington  D.C.;  J.M. Mitchell, 
Jr.  suggested  cooperating  scientist). 

(iii)  Exchange  climatological  data  concerning  the 
energy  balance  of  the  earth-atmosphere  sys- 
tem and  its  geographical  distribution. 

The  USSR  side  would  contribute  surf ace -based 
data  on  the  components  of  heat  balance.  The 
US  side  would  contribute  satellite-based  data 
on  the  planetary  heat  balance,  including 
cloudiness. 
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(iv)  Oceanol^gic  indices  of  scosonal  and  year-to- 
year  variations  of  conditions  in  the  Horth 
Atlantic  Ocean. 

- Measurer,  of  geographical  distribution  of 
ocean-surface  temperature  and  heat  content 
oi  mixed-layer  of  ocean  (by  Institute  of 
Oceanology,  Moscow;  V.G.Kort,  suggested 
cooperating  scientist). 

“ Data  on  oceanic  conditions,  including  month- 
ly average  surface  temperatures  and  their 
variations  over  a period  of  many  years,  for 
^■11  available  locations  m the  ocean  (by 
W°AA  Environmental  Data  Service, * Washington, 
D.C.;  J.Ai.Mj  to  hell,  Jr.,  suggested  coopera- 
ting scientist). 

b)  Millenium  Program  (0  - 1,000  years  before  present) 
This  will  document  in  greater  detail  the  course  of 
global  climate  change  in  earlier  centuries,  thus 
lending  perspective  to  the  course  cf  climate  in  the 
20th  century  in  the  context  of  historical  events 
including  the  "Little  Ice  Age", 

c)  Ice-age  History  Program  (0-p0,000  years  before 
present,  and  older) 

This  will  help  improve  the  understanding  of  the 
oxti'emes  of  climate,  and  of  the  nature  of  the  tran- 
sitions of  climate,  both  within  and  following  a ma- 
jor glacial  event  on  earth;  and  to  clarify  the  geo- 
graphical pattern  of  conditions  and  events  during 
each  phase  of  the  glaciation  and  deglaciation,  as  , . 

needed  to  evaluate  the  accuracy  of  climate  simulation 
models. 

2. 2 Exchange  of. Data 
a)  Tree-ring  widths 

un-noroalized  ring-width  series  f;om  replicated 
trees  in  each  of  11  or  12  subpolar  sites  in  USSR, 
from  the  Finland  border  eastward,  to  the  Kamchatka 
Peninsula  (by  the  Lithuanian  Academy  of  Sciences, 
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Kaunas,  Dr.  Bitvinskas , suggested  cooperating  sci- 
entist; and  by  the  Botanical  Institute,  Leningrad, 
Dr.  Lovelius,  suggested  cooperating  scientist). 

- Ring-width  series  and  detailed  information  concer- 
ning tranofor  functions  relating  ring-width  series 
to  climate,  for  any  or  all  of  the  more  than  40 
available  sites  in  Ilorth  America  as  desired  (by 
University  of  Arizona,  Tucson;  II.C.FrittSj  s gges- 
ted  cooperating  scientist). 

b)  Exchange  pnleoclimatographic  data  acquired  according 
to  the  plans  approved  at  the  Palcoclimatographic  Con- 
ference described  in  2.3.3. 

2.3.3  Action  Proposed. 

a)  Conference  on  pal c o c 1 ima t o g raphyt  a technical  confe- 
rence to  be  held  approximately  June,  1975;  to  be  or- 
ganized to  accomplish  the  specific  objectives  listed 
below,  and  to  include  as  many  of  the  cooperating  sci- 
entists listed  below  as  possible,  in  addition  to 
other  scientists. 


Country 

Institution 

Sunaested  Participants 

USSR 

Institute  of  Geography, 
Moscow 

t*».  I .lleustadt  ,N.S  .Chebotarcva, 
M.G.Grossvald,  V.P.Gruchuk, 

A. A, Velichko,  N.A.Khotinsky 

USSR 

Institute  of  Geology, 
Moscow 

K.N. Nikiforova,  N.V.Kind, 
R.E.Gi German,  E.Korcnova, 
M.A.Pavzner 

USSR 

All-Union  Geological 
Institute  (VSEGEI) 
Leningrad 

I. I. Krasnov 

USSR 

Institute  of  Geology, 
Novosibirsk 

V.N.Saks 

USSR 

Institute  of  Geology, 
Murmansk 

S. A .Strelkov 

US 

Lamont-Doherty  Geologi- 
cal Observatory, 
New-York 

W.S .Broecker,  N.Ondvke 

! 


* To  be  held  in  Russia.  See  page  25. 
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US 

University  of  Wisconsin 

J.Kutzbach 

US 

Brown  University, 
Providence,  Rhode  Island  • 

T.Webb,  J.Imhrie 

US 

University  of  Minnesota 

H. Wright 

US 

I ale  University, 

New  Haven,  Qonnecticut 

M.  Davis 

US 

US  Geological  Survey 

D.Adara,  E. Leopold 

US 

University  of  Maine 

G. Denton 

US 

Ohio  State  University 

T.  Hughes 

US 

University  of  Washington 

S. Port or-  , 

The  objectives  of  the  Conference 

are  as  follows: 

- To  plan  for  the  exchange  of  data  on  C-14-  dated  ice- 
sheet  moraines  from  as  many  locations  as  possible  in 
the  U.S.  and  USSR  (especially  Siberia). 

- To  plan  for  the  exchange  of  historical  palynological 
data  (samples,  pollen  counts,  C-14  dates,  and  j lien 
diagrams)  from  C-1A  dated  stratigraphic  sections 
located  as  widely  as  possible  in  the  U.S.  and  USSR 
including  areas  well  south  of  areas  covered  by  gla- 
cial ice  during  the  Quaternary. 

- To  plan  for  the  exchange  of  modern  pollen  data  (samp- 
les, pollen  counts,  geographic  locations)  from  sites 
distributed  as  widely  as  possible  in  the  U.S.  and  USSR 

- To  exchange  statistical  transfer-function  techniques 
designed,  to  extract  rainfall  and  seasonal  temperature 
estimates  from  fossil  pollen  data, 

- To  assess  the  C-14  dating  facilities  needed  to  accom- 
plish the  first  and  second  objectives  in  the  U.S.  and 

*’  USSR,  and  to  make  appropriate  recommendations  to  in- 
sure the  needed  support. 

- To  formulate  a research  plan  for  acquiring  continuous 
stratigraphic  sections  (both  pollen  diagrams  from  lake 
cores  and  soil  sequences)  covering  all  or  significant 
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portions  of  the  Brunhes  Magnetic  Epoch  (the  pact 
690.000  years)  in  the  U.S.  and  USSR  at  a number  of 
widely-distributed  sites. 

b)  Conference  on  paleoclimatic  node] ling:  a technical 
conference  to  be  held  in  the  U.S.  in  Juno  1976  or 
thereabouts,  designed  to  exchange  results  of  expe- 
riments simulating  selected  past  climates, 

(i)  U.S,  scientists  will  present  results  of  nume- 
rical experiments  aimed  at  simulating  ice-age 
climates  at  the  lust  glacial  maximum  \icing  ge- 
neral circulation  models  (IiCAR,  Boulder,  Colo- 
rado; ^Williams,  W, Washing  ton,  and  R\G.  Barry, 
suggested  cooperating  scientists.  National 
Science  foundation's  lDOK/CLlMAP  Project  Scien- 
tists and  other  institutions  as  appropriate) 

(ii)  USSR  scientists  v.-i.ll  present  results  of  clima- 
te-modelling experiments  carried  on  by  the  In- 
stitute of  Gcograpny,  Moscow  (I.P. Gerasimov, 
suggested  cooperating  scientist)  and  by  the 
Main  Geophysical  Observatory,  Leningrad  (M.I. 
Budyko,  suggested  cooperating  scientist) .In  ad- 
dition, information  on  ice-agc  climate  in  the 
oceans  will  be  presented  by  the  scientific 
staff  of  the  Institute  of  Oceanology,  Moscow 
(A^ysitsin,  suggested  cooperating  scientist) 
and  other  institutions  as  appropriate, 

c)  Exchange  of  scientific  personnel. 

Each  country  shall  have  the  right  to  arrange  for 
working  visits  of  up  to  a total  of  6 man-months 
duration  per  year  by  up  to  5 scientists,  Soviet 
institutions  open  to  visitation  shall  include  tue 
Institutes  of  Geography,  Geology,  and  Oceanology 
in  Moscow,  the  All-Union  Geological  Research  Insti- 
tute in  Leningrad  and  in  Novosibirsk,  and  the  Sibe- 
rian Branch  of  the  USSR  Academy  of  Sciences  in  No- 
vosibirsk, U.S,  institutions  open  to  visitation 
shall  include  those  listed  under  2.3.2.  above. 
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PROJECT  2.  JOINT  STUDIES  OP  THE  EPPECTS  OP  POLLUTION  OP  THE 
ATHOSPHEHE  ON  CLIHATE 

I 

2.4.  Subgroup  A.  Monitor! nr  ntnoanhoTvic  comiti tuents 
and  nnsoBf/ment  of  1 heir  cffcctr?  on 
clir.r.to 

2.4.1.  Scientific  cooperation 

a)  At  present,  the  USA  has  4 baseline  static  ri  and 
10  regional  otationa  with  2 additional  baseline 
otatione  planned  for  the  future.  The  USSR  has  1 
baseline  atatioa  and  5 regional  stations.  It  ic 
desirable  to  enlarge  tho  network  of  stations  and 
to  expand  their  geographical  location.  The  respec- 
tive networks  chould  bo  expanded  us  soon  as  fca— 
Bible  and  largely  completed  by  1S80, 

b)  Both  countries  will  provido  for  tho  accompli- 
ohmont  of  measurements  within  the  WHO  program 
(aoroeol  chemistry  banod  on  the  study  of  the 
atmospheric  turbidity,  rain  chemistry  and  CO^). 
Additional  measurements  beyond  these  minimum 
requirements  are  also  being  made  othoro  are 
planned. 

c)  Both  countries  should  develop  methods  of  mca& 
Burement  and  systems  of  calibration  for  carbon 
dioxide,  ozone,  and  atmospheric  turbidity  (and 
other  optical  properties)  at  tho  baseline  sta- 
tions and  exchange  tho  experience  gained  and 
the  data  from  such  work. 

(i)  The  US  would  provide  a surface  ozone  measu- 
ring system  for  operation  at  a USSR  remote 
Observatory.  The  exchange  would  be  phased 
as  follows: 

1974  - one  electrochemical  concentration 
cell  system 

1975*76  - one  chemiluminescence  (ethylene 
reaction)  system. 

One  precision  ozone  generator. 
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D)  I*  otlr.g  of  Exports  to  Prop are  Dntn  Boso  of  tho  Report 

Clir.ato 


This  group  would  roconnond  how  tho  data  bass  diocuaned  in 
Eeotlon  2^2.2  should  bo  structured,  the  parameters  to  vo  includ- 
ed, their  resolution  in  space  cad  tine,  the  levoi  cf  processing 
of  tho  d..ta  and  the!-'  statistical  treatment.  Tho  group  would 
also  recommend  hen?  such  a data  best  should  bo  compiled  and  on 
appropriate  division  of  responsibility  with  duG  recognition  ox 
tho  nood  for  full  international  cooperation  to  achieve  global 
data  sots  and  tho  rolo  ct  CHRP  in  compiling  global  data  sets 
in  connection  with  GARP  experiments. 

The  group  would  also  recommend  procedures  and  formats  for 
tho  full  exchange  ox  iil^.  data  and  ,?inalyces  that  comprise  the  da- 
ta base  for  tho  present  climate. 

A meeting  of  6 to  8 specialists  '/ill  bo  organised  in  tho 
United  Staton  within  tho  next  12  months.  Dr. Oort  will  orga- 
nise tho  mooting. 

E)  Hooting  o^Rvpojvtn  on  Quaternary^ Glaciation 

This  group  would  study  tho  distributions  and  changes  of  ico 
and  climatic  zones  in  the  USER  during  tho  rocont  uppar-quat9rna- 
Vj  glaciations.  Prof . Inbrio  would  holp  organize  the  meeting 
which  should  be  hold  in  the  USSR  during  19 75. 

3«*.  Long  Duration  Working  Visits  of  Scientists 

It  is  desirable  to  have  scientists  from  each  country  visit 
the  other  country  for  periods  of  6 to  12  months  to  work  on  prob- 
lems connected  with  climatic  variation. 

3»*.1  ‘ The  U.S. would  like  to  issue  invitations  to  Soviet 
•cientlsts  to  visit  the  following  U.S.  institutions  for  6 to  12 
nonth  period. 

(•)  To  cone  to  the  Geophysical  Fluid  Dynamics  Laboratory  of 
IOAA  In  Princeton,  New  Jersey:  (1)  to  study  climate 

sensitivity  and  stability;  (2)  to  work  on  radiation 
effects  of  particulates ; (3)  and  to  study  numerical 

Integration  methods. 
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(b)  To  cone  to  the  National  Contor  for  Atmospheric  Research  • 
in  Boulder,  Colorado  to  study  adjoint  r.othods  of  climatic 
perturbation  and  to  study  effects  of  bob  ico  on  cliaoto. 

(e)  To  come  to  the  University  of  Washington  in  Beattie  to 
study  models  of  ico  dynamics. 

<d)  To  como  to  Stanford  University  to  study  the  analysis 
of  variations  of  temperature  and  pressure  fields  in  the 
uppor  and  lower  stcospbero  during  solar  disturbances. 

The  USSR  would  like  to  issue  invitations  to  U.S.  scientists 
to  visit  the  following  USSR  institutions  for  6 to  12  months  pe- 
riods. 

(a)  To  cone  to  the  Computing  Center  of  tho  Siberian  Branch 
of  the  Academy  of  Sciences  to  study  numerical  modelling 
of  climate. 

(b)  To  come  to  tho  Wain  Geophysical  Observatory  of  tho  Hy- 
dromoteorological  Service  in  Leningrad  to  work  on  prob- 
lems of  measuring  minor  constituents  in  tne  atmosphere. 

(e)  To  come  to  the  Arctic  and  Antarctic  Research  Institute 
in  Leningrad  to  work  on  problems  of  polar  areas  as  they 
effeot  climate. 

(d)  To  come  to  the  Hydrometeorological  Center  in  Moscow 
to  work  ou  problems  of  climate  prediction  and  documen- 
tation. 

For  the  purpose  of  expediting  the  exchange  of  climatic 
data,  a scientist  from  the  USSR  should  visit  the  US 
Rational  Climatic  Center  In  Asheville,  N.C.  At  this 
tine  lists  of  available  climatic  data  and  information 
will  be  exchanged  on  magnetic  tape  formats.  The  Soviet 
scientist  will  bring  an  example  of  a tape  prepared  in 
the  USSR  and  will  be  supplied  with  an  example  of  a tape 
prepared  in  the  U.S.  This  visit  would  probably  last 
one  month  and  mould  bo  carried  out  during  the  next  year. 
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LIST 'OF  PARTICIPANTS 

USA 


" • ’.7.N.HESS 

2.  Dr.S.C.CORONTTI 

3.  Prof . R . CllARLSON 

4.  Dr . E . S . EPS  TSUI 

5.  Mr. 0.0. FLETCHER 

6.  Dr . A . J . GR03SCKER 

7.  Prof . J . 11.31?  IE 

P.  Dr. C.E. LEITH 

„ . i^r.J.M. MITCHELL 

10.  Dr.S.Ll.OLHUICOFF 

11.  Mr. D. II. PACK 

12.  Mrs. F. SETTLE 

13.  Dr. J.SMAGOR  INSKY 

14.  Prof . J .15 . ,7ILC0X 

USSR 


Principal  Delegate,  ROAA/Environmsn 
tal  Research  laboratories 

CIAP/Department  of  Transportation 

University  of  Washington 

NOAA/Environraent  Monitoring  and 
Prediction 

National  Science  Found’ation 

CIAP/Departaent  of  Transpu:  tation 

Department  of  Geological  Sciences, 
Brown  University 

t.ational  Center  for  Atmospheric 
Research 

NOAA/ Environmental  Data  Services 

National  Science  Foundation 

NOAA/Air  Resources  Laooratory 

NOAA/Environment  Research  Labora- 
tories 

NOAA/ Geophysical  Fluid  Dynamics 
Laboratory 

Institute  for  Plasma  Research, 
Stanford  University 


15.  Prof . Ye. P. BORISENKOV 

16.  Prof .M.E.BERLYAND 

•* 

17.  Dr. I. I. BORZENKOVA  • 

18.  Prof  .Li.I.BUDYKO  • 

19.  Dr.A.A.BUZNIXOV 

20.  Prof .0. A.DROZDOV  • 

21.  Prof .L.S.GANDIH 

22.  Dr. A.  A. GRIGORIEV 

23.  Dr. G.G. GROMOVA 

24.  Dr.G.P.GUSHCHIH 


Principal  Delegate,  Voeikov  Main 
Geophysical  Observatory f Leningrad 

Voeikov  Main  Geophysical  Observa- 
tory, Leningrad 

Voeikov  Main  Geophysical  Observa- 
tory , , Leningrad 

Voeikov  Main  Geophysical  Observa- 
tory, Leningrad 

Leningrad  State  University, Leningrad 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad 

Leningrad  State  University .Lenin- 
grad 

Hydrometeorological  Centre  .Moscow 

Voeikov  Main  Geophysical  Observato- 
*7,  Leningrad 
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25.  Dr. I. L. KAROL 

26.  Mr. Yu. E. KAZAKOV 

Dr . G . M . KREKO V 
28  . Dr . V . F .LOGINOV 

29.  Dr.V.P.KELESHKO 

30.  Dr . V . V . MIKHNEVICH 

31.  Prof.E.P.MUSTEL 

32.  Dr.V.V.PENEHKO 

33.  L.R.RAKIPOVA 

34.  Dr.B.I.SAZONOV 

35.  Prof.S.V.SOLONIN 

36.  Dr.V.Ya.SERGIN 

37.  Prof .M.E.SHVETZ 

38.  Dr . V . D . STEPANENKO 

39.  Dr . A . F . TRESHNIKO V 

40.  Dr  .R.F. USMANOV 

41.  Dr.O.B.VASILIEV 

42.  Dr.Z.Ya. VINNIKOV 

43.  Prof  .11. 1. YUDIN 


Voeikov  Main  Geophysical  Observa- 
tory, Leningrad. 

Chief  Administration  of  the  Hydro- 
meteorological Service  of  the  USSR, 
Moscow 

Institute  of  Atmospheric  Optics, 
Tomsk 

All-Union  Research  Institute  of 
Hydrometeoro logical  Information, 
Obninsk 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad 

Institute  of  Applied  Geophysics, 
Moscow 

Hydrometeorological  Centre , Moscow 

Computer  Centre , Novosibirsk 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad 

Voeikov  Main  Geophysical  Observatory 
Leningrad 

Leningrad  Hydrometeorological  Insti- 
tute, Leningrad 

Pacific  Geographic  Institute  of  the 
AS  USSR,  Vladivostok 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad  • 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad 

Arctic  and  Antarctic  Research  Insti- 
tute, Leningrad 

Hydrometeorological  Centre,  Moscow 

Leningrad  State  University, Lenin- 
grad 

Voeikov  Main  Geophysical  Observa- 
tory, Leningrad 

Voeikov  Main  Geophysical  Observato- 
ry, Leningrad 
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CARTOGRAPHIC  ANALYSIS  OF  THE  CHANGING  TIME-SPACE  PATTERNS  IN  THE  HOLOCENE 
POLLEN  RECORD  OF  NORTHEASTERN  NORTH  AMERICA 

Bernabo,  J.  Christopher  and  Thompson  Webb  III,  Dept,  of  Geol.  Sci. , 
Brown  University,  Providence,  R.I.  02912 


Fifty-two  radiocarbon-dated  cores  from  the  northeastern  U.S.  and 
southern  Canada  provided  the  basis  for  mapping  the  distribution  of  sel- 
ected pollen  types  at  each  1,000-year  interval  during  the  Holocene.  The 
sequence  of  these  maps  traces  the  development  of  the  vegetation  over  the 
last  11,000  years.  Three  types  of  maps  were  constructed:  isopollen  maps, 

difference  maps,  and  isochrone  maps.  Isopollen  maps  show  the  distribution 
of  a pollen  type  at  a given  time  and  are  contoured  with  lines  of  equal 
pollen  frequency  (isopolls).  Difference  maps  are  generated  by  subtracting 
successive  isopollen  maps  and  illustrate  the  areas  of  changing  pollen  fre- 
quencies over  selected  time  intervals.  The  production  of  both  of  these 
map. types  converts  a set  of  geographically  disjunct  time-series  (the  pollen 
diagrams)  into  synoptic  regional  patterns.  Recognition  of  these  patterns 
is  crucial  to  deciphering  the  ecological  and  climatic  significance  of  the 
pollen  record.  Isochrone  maps  show  the  "range  limits"  of  specified  polic; 
types  at  1,000-year  intervals  and  are  useful  in  studying  differential 
migration  rates. 

% 

Analysis  of  the  maps  indicates . that  the  major  climatically-induced 
vegetational  changes  of  the  Holocene  culminated  by  8,000  YBP.  Prior  to 
that  time,  changes  in  pollen  distribution  were  closely  related  to  the  posi- 
tion of  the  ice  margin.  By  ca.  8,000  YBP,  the  modem  configuration  of  veg- 
etational formations  was  established  with  conifer-hardwood  forests  in  the 
northern  Great  Lakes,  deciduous  forests  to  the  south  and  prairie  to  the 
west.  Changes  in  the  position  of  the  ecotones  and  in  the  composition  of 
the  pollen  record  were  smaller  in  scale  after  ca.  8,000  YBP  than  before. 
Other  features  of  the  Holocene  pollen  record  depicted  by  the  maps  are  the 
rapid  decline  of  spruce  pollen,  and  westward  extension  of  pine  pollen 
(see  figure)  and  the  mid-Holocene  prairie  expansion  and  retreat. 

• Many  lines  of  independent  evidence  reveal  that  the  early  Holocene  was 
a period  of  rapid  climatic  change.  Comparison  of  the  behavior  of  different 
pollen  types  shows  that  spruce,  pine  and  oak  responded  rapidly  to  these 
changes.  The  shifting  distribution  of  these  genera  appears  to  be  a useful 
monitor  of  the  changing  macroclimate.  Certain  other  genera,  however,  su-h 
as  chestnut,  hickory  and  beech  responded  sluggishly,  migrating  northward 
and  westward  gradually  over  much  of  the  Holocene.  Ecological  factors, 
other  than  climate,  such  as  soil  development,  reproductive  strategy  and 
competitive  interactions  probably  played  a dominant  role  in  controlling 
the  distribution  of  these  genera  (see  M.  B.  Davis,  this  volume). 

The  mapping  of  pollen  data  on  a sub-continental,  scale  emphasizes  the 
broad  geographic  patterns  in  the  vegetation  which  are  indicative  of  macro- 
cllmate  and  changes  in  its  state.  Pollen  maps  also  serve  to  summarize  the 
pollen  record  in  a graphic  form  that  is  readily  comprehensible  to  a wide 
range  of  Interested  Quaternary  scientists,  outside  palynology. 
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PALYNOLOGICAL  EVIDENCE  OF  ECOTONAL  MIGRATIONS  IN  EASTERN  NORTH  AMERICA 

Webb,  Thompson  III  and  JL  Christopher  Bernabo,  Dept,  of  Geo.  Sci., 

Brown  University,  Providence,  R.I.  02912 

Pollen  data  monitor  the  vegetation  through  time  and  space.  Comparison 
of  maps  of  the  modern  pollen  data  with  maps  of  the  modern  vegetation  show 
the  patterns  of  formations  and  forest  types  to  be  represented  clearly  with- 
in the  pollen  data  (Webb,  1974,  Ecology) . Transitions  (or  ecotoner.)  between 
forest  types  or  formations  are  therefore  evident  upon  the  maps  of  modern 
pollen  and  are  recorded  by  changes  in  the  relative  abundances  of  certain 
pollen  types  lrom  one  region  to  the  next.  The  relative  weightings  of  the 
ecological  factors  controlling  the  positions  of  the  various  ecotones  is  as 
yet  unspecified.  Study  of  the  behavior  of  these  ecotones  through  time,  how- 
ever, may  help  in  defining  the  relative  importance  of  the  various  ecological 
factors. 

Bernabo  ejt  al.  (see  Abstract,  this  volume)  have  prepared  maps  showing 
the  distribution  of  pollen  types  at  selected  times  during  the  Holocene  in 
eastern  North  America.  The  sequence  of  these  maps  from  11,000  YBP  to  pre- 
sent shows  that  ecotones  change  in  character  and  position  through  time.  On 
the  modern  map'-,  of  pollen  distribution,  the  ecotone  between  deciduous  for- 
ests and  conifer-hardwood  forests  in  the  Midwest  (which  coincides  with  the 
"tension  zone"  there)  is  marked  by  the  south  to  north  change  from  the  domi- 
nance of  oak  pollen  to  the  dominaned  of  pine  pollen.  This  same  change  in 
pollen  dominance  can  be  used  to  trace  the  position  of  this  ecotone  in  the 
past  (see  figure).  Early  in  the  Holocene,  the  position  of  this  ecotone 
changes  drastically,  probably  in  direct  response  to  the  fast  changing  cli- 
matic patterns  of  that  period.  After  8,000  YBP,  the  position  of  this  eco- 
tone remains  relatively  stable  but  the  composition  of  the  pollen  spectra 
continue  changing  on  both  sides  of  the  boundary  as  new  pollen  types  appear. 
This  relative  stability  in  the  location  of  the  ecotone  is  in  direct  contrast 
to  the  changing  composition  of  the  formations  it  separates,  and  this  stabil- 
ity reflects  a different  weighting  of  ecological  factors  that  control  the 
ecotone' s position  from  the  weighting  of  the  factors  controlling  community 
composition.  Probably  the  former  is  mainly  controlled  by  macroclimatic 
factors,  e.g. , air  mass  patterns,  whereas  the  latter  is  controlled  largely 
by  edaphic,  biological  and  small-scale  climatic  factors. 
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Isoc&rones  of  a pollen-defined  ecotone  in  10s  years.  Present  positon  along 

^-7  contour.  , ■ 
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ON  THE  USE  OF  TREE  RINGS,  POLLEN  AND  MARINE  PLANKTON  IN  RECONSTRUCTING  PAST 
CLIMATES 

Webb,  Thompson  III  and  Douglas  R.  Clark,  Dept,  of  Geol.  Sci.  , Brown 

University,  Providence,  R.I.  02912 

Large  collections  of  tree  ring,  pollen  and  marine  plankton  data  enable 
the  study  of  eco-climatic  interactions  over  wide  geographic  areas,  and 
over  long  time  spans.  Multivariate  statistical  tcchinques  like  those  used 
by  Iritts  (Abstract,  this  volume)  provide  a means  for  obtaining  quantitative 
estimates  of  past  cli.mates  from  the  biological  data.  When  enough  cores 
containing  one  type  of  these  data  are  available  in  a region,  a map  depict- 
ing the  estimated  climate  can  be  prepared  for  selected  time-intervals  with- 
in the  age  range  of  the  cores  and  such  paleoclimatic  maps  are  new  becoming 
available.  6 


.One  remarkable  feature  of  these  three  types  of  data  is  that  they  com- 
plement each  other  in  their  time  resolution  and  their  length  of  record. 

In  the  data  sets  suite ble  for  mapping,  the  annually  spaced  tree  rings  give 
decade  by  decade  summaries  of  the  past  500  years;  pollen  data  that  are 
generally  sampled  in  200-year  intervals  give  1,000  year  by  1,000  year  sum- 
maries of  the  past  12,000  years;  and  marine  plankton  data  sampled  in  ca. 
3,000-year  intervals  give  10,000  year  by  10,000  year  summaries  of  the  past 

150,000  years.  , 

• ' • 

Overlap  between  the  three  different  data  sets  also  occurs.  For  in- 
stance, pollen  from  annually  laminated  lake-sediments  can  have  the  time 
precision  of  tree  rings,  and  tree  ring  sequences  can  be  extended  the 
length  of  Holocene  pollen  records.  This  overlap  is  important  because  1) 
where  two  data  sets  occur  in  one  locality,  climatic  estimates  derived  from 
one  data  set  can  be  used  to  verify  the  estimates  from  another  data  set,  and 
2)  where  the  two  data  sets  occur  in  different  localities,  possible  telecon- 
ncc.tions  in  the  climates  of  the  two  regions  can  be  studied.  An  opportunity 
of  the  latter  type  exists  given  the  mapped  pollen  records  of  eastern  North 
America  (Bernabo  j^t  al.  , Abstract,  this  volune)  and  the  mapped  foraminifera 
data  of  Ruddiraan  (Abstract,  this  volume)  in  the  North  Atlantic.  These  two 
data  sets  provide  a basis  for  studying  the  interaction  between  climatic 
patterns  over  land  and  sea  during  the  early  Holocene. 

In  eastern  North  America  between  10,000  and  9,000  B.P.,  the  "late- 
glacial  boreal  forest"  largely  disappears  and  high  values  of  pine  pollen 
appear  within  100  km  of  the  Laurentide  ice  margin.  This  change  indicates 
a steepening  of  the  climatic  gradient  south  of  the  ice  margin  and  creates 
a meteorological  situation  conducive  to  strong  westerly  flow  near  the  ice 
margin.  This  proposed  circulation  pattern  supports  the  one  hypothesized 
by  Ruddiman  to  explain  the  distribution  of  sea-surface  temperatures  in  the 
North  Atlantic  at  9,300  B.P. 
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PALEOECO LOGICAL  SIGNIFICANCE  OF  THE  CHANGING  PATTERNS  IN 
THE  HOLOCENE  POLLEN  RECORD  OF  NORTHEASTERN  NORTH  AMERICA 
Bernnbo,  J.  Christopher  and  Webb,  Thompson,  III,  Do- 
par  I'.incnL  of  Geological  Sciences,  Brown  University, 
Providence,  R.i.  02912 

■laps  showing  the  distribution  of  pollen  from  the  major  tree 
genera  and  herbs,  for  each  1,000  yr.  interval  during  the 
lolocenc,  were  constructed  using  52  radiocarbon-dated 
cores.  The  production  of  such  maps  converts  a set  of  geo- 
graphically disjunct  time  series  (the  pollen  diagrams)  into 
synoptic  regional  patterns.  Examining  the  dynamics  of 
these  patterns  is  crucial  to  interpreting  the  climatic  and 
ecologica.l  significance  of  the  pollen  record.  Analysis  of 
the  maps  indicates  that  the  major  climatically-induced  veg- 
etational  displacements  of  the  Holocene  occurred  before 
3, 000  YBP.  Prior  to  that  time,  changes  in  pollen  distribu- 
tion were  closely  related  to  the  position  of  the  shifting 
ice  margin.  By  8,000  YBP  the  modern  configuration  of  vcgc- 
tational  regions  was  established  with  mixed  conifer— hard- 
wood forests  in  the  northern  Great  Lakes,  deciduous  forests 
to  the  south  and  prairie  to  the  west.  Changes  after  8,000 
YBP  were  smaller  in  scale  and  altered  the  areal  extent  not 
the  gross  configuration  of  the  vegctational  regions.  Com- 
parison of  the  behcivior  of  different  genera  reveals  that 
spruce,  pine  and  oak  responded  rapidly  to  changing  environ- 
mental conditions  and  arc  useful  paleoclimatic  indicators. 
However , hickory,  hemlock,  beech  and  chestnut  gradually 
migrated  northward  during  much  of  the  Holocene.  The  slower 
climatic  responses  of  these  genera  makes  their  paleoenvir- 
onmental  message  unclear. 
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The  relationship  between  oceanic  and  continental  paleo— 
climates  is  unclear;  adequate  data  on  air/sea  interaction  being 
available  only  for  the  past  several  hundred  years.  Fine  scale 
correlation  of  marine  and  terrestrial  paleoclimatic  data  is 
limited  by  the  often  brief  and  abrupt  nature  of  climatic  change 
and  by  errors  associated  with  age-dating  techniques.  In  order 
to  overcome  these  correlation  problems,  we  analyzed  marine  and 
terrestrial  climatic  indicators  from  identical  samples  in  C^-4 
dated  marine  cores  thereby  obtaining  a direct  correlation  of 
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climatic  events. 

Both  pollen  and  planktonic  foraminifera  were  analyzed  from 
cores  V24-1  and  V26-176  taken  from  the  continental  rise  off  the 
mouth  of  Chesapeake  Bay  (Figure  1)  . The  sediment  in  these  cores 
is  an  olive-gray  lutite  with  coarse  fraction  (weight  greater  than 
63  microns)  rarely  exceeding  10  percent.  The  carbonate  content 
of  the  cores  fluctuates  from  4 to  32  percent  with  higher  values 
being  characteristic  of  the  last  8000  years  (Figure  2) . In 
these  cores  C-^4  dates  were  used  in  conjunction  with  carbonate 
minima  and  maxima  and  floral  datum  planes  to  establish  a down- 

core  chronology  and  correlation  between  cores. 

% 

Pollen  concentration  varied  from  800  to  12,000  grains/gram 
thereby  permitting  identification  of  a significant  number  of 
pollen  grains  in  each  sample.  Pollen  diagrams  from  these  marine 
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cores  exhibit  the  following  features  (Figures  3 and  A) . Boreal 
vegetation  was  replaced  by  thermophilous  forest  elements  about 
10,000  YBP,  the  late  glacial/Holocene  boundary.  Glacial  (before 
10,000  YBP)  pollen  assemblages  are  characterized  by  pine,  spruce, 
fir,  and  herbs,  which  in  the  Early  Holocene  are  supplanted  by  a 
pine-oak  assemblage  in  which  hemlock  is  briefly  abundant.  Pol- 
len from  deciduous  trees  continues  to  increase  in  abundance  until 
3500  YBP  when  oak  declines  and  pine  re-expands.  Comparison  of 
these  marine  pollen  diagrams  to  pollen  diagrams  from  the  Chesa- 
peake drainage  basin  demonstrates  that  pollen  assemblages  from 
these  marine  cores  accurately  reflect  the  timing  of  vegetation 
changes  on  land. 

Planktonic  foraminifera  were  analyzed  using  transfer  function 
F13  to  obtain  down-core  estimates  of  paleo-oceanographic  para- 
meters. In  addition,  because  these  cores  are  located  in  waters 
influenced  by  the  Gulf  Stream,  certain  species  of  planktonic  for- 
aminifera were  used  as  Gulf  Stream  indicators.  In  particular, 

10  percent,  or  greater  abundance  of  Globogernoides  sacculifer 
s appears  to  be  restricted  to  the  average  position  of  the  warmest 
part  of  the  Gulf  Stream.  Down  core  estimates  of  winter  tempera- 
ture indicates  that  the  maximum  temperature  .in  this  area  was 
reached  about  8000  YBP  (Figure  5) . Furthermore,  the  increased 

abundance  of  G.  _sacculifer  suggests  that  this  temperature  maxi- 
> « 
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mum  is  associated  with  the  movement  of  the  Gulf  Stream.  From 
8000  YBP  to  the  present,  both  cores  exhibit  generally  decreasing 
surface  temperatures  (Figure  5)  related  to  the  increased  import- 
ance of  the  subpolar  fauna  (Figure  6) . Today  the  dominant  faunal 
assemblage  at  the  site  of  both  cores  is  subpolar  (Figure  7) 
However,  V26-176  is  very  close  to  the  boundary  between  the  sub- 
polar and  tropical/gyre  margin  assemblage.  Before  8000  YBP,  the 
two  cores  show  a diachronous  pattern  of  estimated  winter  temper- 
ature changes  related  to  the  complex  interaction  of  the  subpolar 
and  tropical  fauna  (Figure  6) . 

Comparison  of  foraminifera-derived  temperature  estimates 
and  poilen  fluctuations  in  these  cores  suggests  that  in  the  area 
studied,  terrestrial  climatic  events  lag  somewhat  behind  changes 
in  the  ocean  (Figure  8) . However,  terrestrial  changes  appear  to 
be  more  closely  related  to  variations  in  slope  water  temperature 
than  to  Gulf  Stream  migrations. 


- D4  - 

Figure  Captions 

Figure  1.  Sketch  map  showing  location  of  cores  and  study  area. 

The  locations  and  depths  of  the  cores  are  as  follows: 

V24-1,  36°30'N  lat. , 73°30'W  long.,  3012  m. ; V26-176,  36°03'N 
lat.,  72°23'W  long.,  3542  m. 

Figure  2.  Down  core  variations  in  calcium  carbonate  and  cor- 
. relation  points  indicated  by  C (carbonate  correlation) , D 
(radiocarbon  date)  , and  F (floral  datum) . Tne  correlation 
points  are  as  follows:  C^,  carbonate  maxima  C^-4-dated  at 

about  8000  YBP  in  both  cores;  C2,  carbonate  minima  co- 
occurring with  F2;  C3,  first  carbonate  maxima  below  C2 ; Dj , 
8050  + 250  C14  YBP;  D2 , 12,400  + 350  C14  YBP;  D3,  3830  +.140 
C14  YBP;  D4,  8180  + 200  C*4  YBP;  F1#  floral  datum  1;  F2, 
floral  datum  2;  F3,  floral  datum  3 (floral  datum  planes 
identified  on  Figures  3 and  4) . Letter  with  identical  sub- 
scripts are  judged  to  be  the  same  age  and  are  used  to 
establish  time  planes  between  the  cores. 

Figure  3.  Pollen  diagram  from  core  V24-1  showing  relative  fre- 
quency of  selected  pollen  taxa.  Note  that  the  percent 
scale  for  pollen  taxa  changes;  Pinus,  Quercus,  and  Picea 


are  at  one  scale  and  the  remaining. taxa  at  another. 
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Figure  4.  Pollen  diagram  from  core  V26-176  showing  relative  | 

frequency  of  selected  pollen  taxa . Variations  in  percent 
* 

scale  are  the  same  as  for  Figure  3. 


Figure  5.  Down  core  estimates  of  winter  sea-surface  temperature 
(TW)  and  percent  abundance  Globogerinoides  sacculifer . 

Figure  6.  Down  core  variations  in  selected  planktonic  foramini- 
feral  assemblages.  Factor  loading  value  indicates  the  re- 
lative importance  of  each  assemblage  in  a sample. 


Figure  7.  Present  distribution  of  dominant  assemblages  off  Cape 
Hatteras,  North  Carolina. 

Figure  8.  Direct  correlation  of  land/sea  climatic  changes. 
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CHANGES  IN  CLIMATIC  REGIME 
NEAR  CAPE  HATTERAS,  N.C. 
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ABSTRACT 

Transfer  functions  are  empirical  equations  for  making 
quantitative  estimates  of  past  atmospheric  and  oceanic  con- 
ditions from  paleontologic  data.  They  are  calibrated  on  a 
set  of  geographically-distributed  observations  which  record 
both  the  modern  climate  and  the  biotic  response  to  climate 
as  reflected  in  paleontological  samples.  The  technique  has 
been  used  on  tree  rings  and  pollen  as  well  as  on  planktonic 
microfossils  (including  foraminif era,  coccoliths,  and  radio- 
laria) . Down-core  applications  in  deep-sea  Quaternary  sedi- 
ments yield  time  series  reflecting  local  changes  in  seasonal 


temperatures  and  salinity.  When  applied  to  samples  from  a 
chronostratigraphic  horizon,  transfer  functions  yield  synop- 
tic charts  of  past  climates.  In  this  form  paleontological 
data  are  readily  interpretable  by  meteorologists  and  physi- 
cal oceanographers , and  can  be  entered  as  boundary  conditions 
into  numerical  models  of  climate. 
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INTRODUCTION 


Since  the  early  years  of  the  last  century,  when  Esmark 
and  Agassiz  made  their  astonishing  discovery  of  ancient  ice 
ages,  the  story  of  past  climates  has  been  elaborated  in  con- 
siderable detail.  This  story  has  been  documented  from  a 
wide  variety  of  sources,  ranging  from  moraines  which  fix 
the  position  of  ancient  ice  bodies  to  studies  of  fossil  pol- 
len and  plankton  which  provide  panoramas  of  past  continents 
and  oceans.  Yet  these  findings  remain  unexplained  in  terms 
of  verified  theories  of  climatic  change.  Why? 

Several  basic  reasons  may  be  cited,  not  least  among 
them  the  extraordinarily  complex  structure  of  the  climate 
system.  The  heart  of  the  problem,  however,  is  the  diffi- 
culty of  achieving  an  effective  interplay  between  climatic 
theory  and  paleoclimatic  data.  This  difficulty  results 
from  a mismatch  in  the  time-space  structure  of  the  observa- 
tional and  theoretical  realms.  Where  dynamic  theory  is 
strong  i.e.,  in  its  ability  to  explain  synoptic  patterns 

of  climate  in  terms  of  atmospheric  and  oceanic  circulation  

the  corresponding  sets  of  paleoclimatic  data  on  spatial  pat- 
terns of  climate  are  sparse.  And  where  kinetic  theory  is 
weak  — — i.e. , in  its  ability  to  analyze  climatic  change  — — 
we  have  an  abundance  of  geologic  time  series  documenting 
climatic  change.  This  mismatch  is  now  being  overcome. 
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firstly  by  translating  geological  observations  into  quanti- 
tative estimates  of  those  climatic  variables  which  physical 
theories  use  and  predict;  and  secondly,  by  producing  synop- 
tic maps  of  the  paleoclimatic  estimates  in  areas  where  sam- 
ples are  densely  arrayed  and  adequately  dated. 

Transfer-function  techniques,  which  are  the  subject  of 
this  paper,  are  playing  a major  role  in  this  task.  These 
techniques  use  multivariate  statistical  procedures  to  ex- 
tract quantitative  estimates  of  air  and  sea-surface  proper- 
ties from  the  fossil  records.  If  the  stratigraphic  coverage 
permits,  these  quantitative  estimates  can  be  mapped  for  sel- 
ected time  intervals  during  the  Quaternary  (Fritts  et:  al_. , 
1971;  McIntyre  eit  aj^. , 1975).  In  this  form  the  data  can 
then  be  used  either  as  boundary  condition  for  running  mathe- 
matical models  of  the  global  circulation;  or  as  a basis  for 
verifying  the  results  of  such  models . Although  still  in 
their  infancy,  these  powerful  models  are  the  main  focus  of 
a growing  scientific  effort  to  understand  past,  present, 
and  future  climates. 

This  paper  summarizes  some  of  the  proceedings  of  a con- 
ference on  transfer  functions  held  at  the  University  of  Wis- 
consin, U.S.A.,  April  3-5,  1974,  and  funded  by  special 
grants  from  the  ARPA  Climate  Dynamics  Program  and  the 


[ 

I 


r 


i 


[ | 

II  1 


j 

* 


National  Science  Foundation's  International  Decade  of  Ocean 
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Exploration.  The  initiative  for  tbe  Wisconsin  conference 
stemmed  from  a 1973  Report  of  SCOR  Working  Group  40,  which 
recommended  a conference  that  would  review  the  various 
models,  analyze  their  compatibility,  assess  their  strengths 
and  weaknesses,  evaluate  the  need  for  further  methodologi- 
cal work,  and  produce  reports  to  promote  wider  use  and 
understanding  of  these  techniques. 

The  objective  of  this  paper  is  to  make  available  to 
the  community  of  scientists  studying  the  oceanic  micro- 
paleontological  record  an  overview  of  the  subject  of  trans- 
fer functions.  We  give  references  to  previous  work;  des- 
cribe briefly  the  model  used  in  deriving  transfer  functions 
from  micropal eon to logic  data;  and  outline  an  algebraically 
homogeneous  set  of  methods  which  have  been  used  to  apply  to 
the  model.  In  another  paper  (Webb  ejt  al^. , in  preparation) 
we  flesh  cut  the  algebraic  details  of  these  methods,  pre- 
sent the  results  of  a series  of  algebraic  experiments  de- 
signed to  explore  their  properties  and  make  specific  recom- 
mendations for  research  to  solve  some  basic  problems. 
PREVIOUS  WORK 

Transfer  function  techniques  were  first  applied  to 
marine  plankton  data  by  Imbrie  and  Kipp  (1971)  , to  terres- 
trial pollen  data  by  Webb  and  Bryson  (1972)  , and  to  tree 
ring  data  by  Fritts  et  al . , (1971).  Recent  applications 
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and  improvements  of  these  techniques  of  special  interest 
to  marine  micropaleontologists  appear  in  articles  by 
Gardner  and  Hays  (1975),  Hecht  and  Kipp  (1974),  Imbrie  et 
al.,(l973>,  Kipp  (1975),  Lozano  and  Hays,  (1975),  McIntyre 
et.  al. , (1975),  Prell  et  al.,  (1975),  Sachs  (1973),  and 
Sancetta  et  al . , (1973). 

THE  BASIC  MODEL 
Definitions 

Let  the  matrix  X be  a defined  set  of  biological  re- 
sponse properties,  e.g.,  relative  abundances  of  organisms 
measured  over  a defined  realm  of  space  and  time  (x,y;t). 
Let  C be  a set  of  physical  variables  of  climate,  either 
marine  or  atmospheric,  measured  over  the  same  time-space 
realm  and  assumed  to  be  causally  related  to  X.  Examples 
include  seasonal  sea-surface  and  air  temperatures  and  sea- 
sonal values  of  precipitation  and  salinity.  Let  D be  a 
set  of  other  physical  variables  of  the  system  which  to- 
gether with  C completely  determine  X.  D often  includes 
such  dif f icult-to-measure  factors  as  nutrient-availability, 
soil— texture , dissolution,  and  anthropogenic  disturbance. 
Thus  the  model  is  completely  deterministic,  in  the  sense 
that  the  observations  in  X are  explained  as  responses  to 
the  ambient  field  of  physical  variables.  No  allowance  is 
made  for  the  operation  of  random  factors,  or  of  historical 
influences  (relict  biotas,  migration  patterns,  etc.). 
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Then,  if  D = O or  is  constant,  the  system  consists  of 
X,  C,  and  a set  of  ecological  response  functions  Re: 

(1)  X = Re  (C)  . 

If  D is  nonconstant,  we  must  consider  the  total  response 
function  Rt: 

(2)  X = Rt  (C:D)  . 

A fundamental  problem  of  paleoclimatology  is  to  find  a set 
•of  transfer  functions  0 such  that  C can  be  estimated  given 
X,  as  follows: 

(3)  C = 0 (X) 

In  practice  0 is  obtained  by  direct  empirical  methods  and 
not  by  inversion  of  Re  or  R^.  Although  such  inversion  pro- 
cedures would  have  many  theoretical  advantages,  and  should 
be  explored  further,  research  to  date  has  uncovered  many 
practical  problems. 

The  class  of  empirical  solutions  to  equation  (3)  that 
will  be  discussed  here  take  the  form 

(4)  C = XB 

in  which  B is  a matrix  of  empirical  calibration  functions . 
Thus  B represents  one  subclass  of  0,  and  relates  to  a moni- 
toring system  with  domain  X and  range  C.  The  X and  C used 
to  derive  the  calibration  functions  B are  the  calibration 
data-set  (Xm»Cm) . The  X to  which  the  calibration  functions 
are  applied  is  the  application  data  set  (Xf) . Xm  is  a spa- 
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tial  array  of  observation  points  over  the  modern  land  sur- 
face or  seabed;  Xf  is  a temporal  array  of  fossil  data,  e.g. 
micropaleontological  data  from  a core. 

As  defined  in  Imbrie  et  al. , (1973,  p.  11),  equations 
that  relate  biological  indices  to  various  properties  of  the 
ocean  or  atmosphere  "may  be  termed  transfer  functions  in 
the  sense  that  they  are  a means  of  processing  one  time- 
varying  signal  (or  set  of  signals)  in  a core,  to  yield 
another  signal  consisting  of  paleotemperature  estimates." 

All  of  the  paleontological  transfer  functions  written  to 
date  either  are,  or  can  be  reduced  to,  simple  linear  trans- 
formations symbolized  as  B in  equation  (4)  . To  emphasize 
this  fact,  and  to  encourage  the  application  of  more  soph- 
isticated procedures,  we  have  chosen  to  define  B as  a set 
of  calibration  functions  and  to  consider  B as  a subset  of 
transfer  functions  (0)  . 

If  it  were  possible  to  derive  0 deductively  from  a 
process-oriented  model  based  on  ecological  first  princi- 
pals, the  predictive  power  of  this  model  would  be  free  of 
many  of  the  problems  encountered  in  its  use  (see  Webb  et 
iJL*  * i-n  prep.).  Given  the  current  state  of  ecological 
theory,  however,  the  empirical  approach  represents  the  only 
practical  solution  to  equation  3,  and  we  are,  therefore, 
stuck  with  its  limitations  and  have  no  choice  but  to  con— 


Imbrie  and  Webb 


E8 


front  the  problems  arising  from  its  application. 

Major  Assumptions 

When  0 is  applied  to  data  on  fossils  several  assump- 
tions must  be  made,  as  outlined  in  Imbrie  and  Kipp  (1971, 
p.  79)  and  in  Webb  and  Bryson  (1972,  p.  74) . The  four  most 
important  are: 

1)  That  the  ecosystem  under  study  has  not  changed 
significantly  during  the  interval  represented  by  the  appli- 
cation data  set  Xf . Specifically,  the  assumption  is  made 
that  the  species  recorded  in  Xm  are  essentially  the  same 
biological  entities  as  those  in  Xf,  and  that  they  have  not 
changed  significantly  their  ecological  responses  to  indivi- 
dual physical  properties  of  the  environment.  This  assump- 
tion is  a highly  restricted  form  of  the  principle  that  the 
present  is  the  key  to  the  past  and  requires  that  contempo- 
rary spatial  patterns  can  be  used  to  interpret  changes 
through  time. 

2)  That  the  biota  represented  by  X are  systematically 
related  to  the  physical  attributes  of  the  medium  in  which 
they  dwell. 

3)  That  the  variables  recorded  in  C are,  or  are 
linearly  related  to,  ecologically  significant  aspects  of 
the  ocean  or  the  atmosphere. 

4)  That  mathematical  equations  representing  linear 


Imbrie  and  Webb 


E 


combinations  of  biological  taxa  can  reflect  the  biotic 
responses  to  physical  changes  adequately  enough  to  yield 
accurate  calibration  functions. 

ALGEBRAIC  METHODS 

Several  multivariate  methods  can  be  used  in  the  cal- 
culation of  the  calibration  functions.  These  techniques 
include  multiple  regression,  stepwise  multiple  regression 
principal  components  analysis  plus  multiple  regression 
(Imbrie  and  Kipp,  1971) , canonical  correlation  analysis 
(Webb  and  Bryson,  1972) , and  distance  coefficients  plus 
regression  (Hecht,  1973) . These  techniques  are  described 
and  their  results  are  compared  in  Webb  et  al.  (in  prep.) 
and  in  Hecht  and  Kipp  (1974) . 
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APPLICATIONS 

Transfer  functions  were  first  applied  to  planktonic 
f oraminif era  in  a deep-sea  Caribbean  core  (Imbrie  and 
Kipp,  1971) . Since  then  the  technique  has  been  used  on  a 
number  of  planktonic  taxa  in  Quaternary  cores  raised  from 
the  floors  of  many  oceans.  Down-core  applications  include 
studies  of  radiolaria  from  the  North  Pacific  (Sachs,  1973a 
and  1973b),  Eastern  Pacific  (Moore,  1973),  and  Antarctic 
(Hays  et  al . , 1975;  Lozano  and  Hays,  1975) ; studies  of 
planktonic  foraminifera  from  the  Norwegian  Sea  (Kellogg, 
1975),  North  Atlantic  (Sancetta  jrt  al^  , 1973),  and  Carib- 
bean (Prell  and  Hays,  1975;  Prell  eb  a_l . , 1975),  Equatorial 
Atlantic  (Gardner  and  Hays,  1975),  Gulf  of  Mexico  (Brunner 
and  Cooley,  1975),  and  Pacific  (Luz,  1973);  and  studies  of 
coccoliths  from  the  North  Atlantic  (Roche  et  a_l.  , 1975)  , 
and  Pacific  (Geitzenauer  et  al . , 1975).  In  addition, 
work  in  progress  by  N.  Maynard,  L.  Burekle,  C.  Sancetta, 
and  H.  Schrader  prove  that  Quaternary  diatoms  can  be 
interpreted  by  transfer  function  techniques. 

Most  of  the  studies  just  cited  give  transfer  func- 
tions for  estimating  three  different  physical  properties 
of  the  surface  water:  summer  temperature,  winter  temper- 

ature, and  salinity.  Although  these  properties  are  highly 


Imbrie  and  Webb 


E 11 


correlated  within  modern  ocean  basins,  downcore  results 
show  considerable  independence.  VJhereas  thermal  and 
salinity  estimates  show  a high  positive  correlation  at 
mid-  and  high-latitude  sites  in  the  North  Atlantic  during 
the  past  130,000  years  (Figure  1) , a negative  correlation 
exists  between  surface-water  temperatures  and  salinities 
estimated  at  many  low  latitude  sites  (Figure  2)  . The 
positive  correlation  in  the  higher  latitudes  is  attributed 
to  changes  in  circulation  pattern  (McIntyre  et  a]_  , 1972)  : 
cold,  low-salinity  water  penetrates  southward  during  gla- 
cial phases  and  replaces  warm,  higher-salinity  water. 

The  decrease  in  salinity  associated  with  this  southward 
penetration  of  polar  waters  apparently  exceeds  the  increase 
that,  for  the  average  ocean,  must  result  from  the  growth 
of  ice  sheets.  The  negative  correlation  between  tempera- 
ture and  salinity  estimates  observed  at  many  low-latitude 
sites  (Figure  2)  must  be  due  to  some  combination  of  two 
effects  which  would  act  in  the  direction  of  increasing 
salinity  during  a glacial  phase:  the  global  build-up  of 

continental  ice,  and  local  increases  in  the  evaporation- 
precipitation  ratio  (Prell  _et  ^1 . , 1975). 

How  accurate  are  paleoclimatic  estimates  derived  from 
transfer  functions?  The  first  step  in  answering  this 
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fundamental  question  is  to  evaluate  bow  transfer  funccions 
perform  on  various  sets  of  data  reflecting  modern  conditions, 
After  calculating  the  standard  error  of  estimate  in  the  cal- 
ibration data  set  (Figure  3; , the  performance  of  the  trans- 
fer function  can  be  tested  on  modern  samples  representing 
different  geographic  areas.  Kipp  (1975),  for  example,  has 
carried  out  such  a test  for  a transfer  function  calibrated 
on  a North  Atlantic  data  set.  When  applied  to  60  samples 
from  the  South  Atlantic  sea  bed,  the  estimates  are  accurate 
within  the  previously  calculated  confidence  interval.  Al- 
though such  a result  does  not  in  itself  validate  down-core 
estimates,  it  does  encourage  belief  in  estimates  derived 
from  late-Qua ternary  Atlantic  sediments. 

The  second  step  in  verifying  down-core  estimates  is 
to  screen  out  any  sample  whose  taxonomic  composition  or 
dissolution  level  lies  outside  the  range  represented  in 
the  calibration  data  set.  Such  samples  represent  condi- 
tions for  which  no  modern  analog  is  known. 

The  third  step  in  the  procedure  for  verifying  down- 
core  estimates  is  the  comparison  of  the  results  of  any 
transfer  function  with  independent  paleoclimatic  estimates. 
Concordant  estirites  encourage  belief;  discordant  estimates 
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discourage  belief,  but  may  raise  important  quesrions  for 
further  research. 

One  class  of  independent  estimates  results  from 
the  application  of  transfer  functions  to  different  bio- 
tic groups.  Another  class  of  independent  estimates  is 
based  on  isotopic  methods.  Although  applications  of  the 
oxygen-isotope  pal  eotemperature  method  must  now  take  into 
account  the  discovery  by  Shackleton  and  Opdyke  (1973) 
that  the  water-composition  effect  dominates  the  thermal 
effect,  these  authors  outline  an  alternate  procedure  for 
obtaining  temperature  estimates.  When  this  alternati  - 
procedure  wes  applied  to  one  Caribbean  core  (Shackleton 
and  Opdyke,  19/3),  the  result  was  concordant  with  the 
previous  estimate  made  by  transfer  function  methods  (Imbrie 

— — * ' 1973).  Clearly,  much  more  work  of  this  kind  re- 
mains to  be  done. 

After  completing  the  testing  of  the  estimates  just 
described,  the  transfer  function  techniques  can  be  used 
in  constructing  paleoclimatic  maps.  In  this  form 
paleoenvironmental  data  are  readily  interpretable 
by  meteorologists  and  physical  oceanographers,  and  as 
discussed  above  can  be  entered  as  boundary  conditions  into 


Imbr' 3 and  Webb 


E 14 


numerical  models  of  climate.  One  goal  of  the  International 
Decade  of  Ocean  Exploration's  CLIMAP  Project  is  to  prepare 
temperature  and  salinity  charts  for  past  oceans  and  arrange 
for  their  analysis  by  numerical  models.  An  example  of  such 
a map  is  given  in  Figure  4,  where  the  sea -surface  tempera- 
ture for  the  North  Atlantic  Ocean  during  an  ice-age  summer 
is  reconstructed. 

Perhaps  the  most  important  feature  of  this  ice-age  re- 
construction is  the  demonstration  of  a structured,  synoptic 
pattern  of  climatic  change.  Portions  of  the  subtropical 
gyre,  for  example,  were  actually  slightly  warmer  during  the 
glacial  phase  than  they  are  today.  Elsewhere  the  anomaly 
is  negative;  but  while  the  Caribbean  cools  only  1°C  to  2°c, 
the  temperature  change  in  many  higher  latitudes  reaches  and 
in  some  places  exceeds  10°C.  By  themselves,  however,  neither 
the  small  changes  in  the  low  latitudes  nor  the  large  changes 
in  mid-latitudes  provide  an  adequate  record  of  the  climatic 
change.  Any  dynamic  explanation  of  the  change  in  the  North 
Atlantic  must  address  uhe  entire  synoptic  pattern.  It  is 
here  that  transfer  functions  can  make  their  greatest  con- 
tribution to  the  study  of  paleoclimates : normalizing  to  a 
common  temperature  scale  paleoenvironmental  inferences 
based  on  different  groups  of  fossils. 
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Figure  Captions 

Figure  1.  Faunal  indices  Ts,  Tw,  and  S vs  depth  in 
North  Atlantic  core  V23-82.  Thick  bars  indicate  80%  con- 
fidence intervals  for  analytical  precision  of  indices. 
Horizontal  lines  indicate  samples  lacking  coccoliths 
(McIntyre  <?t  a_l . , 1972)  . Pattern  in  column  indicates  sam- 
ples with  high  mineral  detritus  content  in  the  > 149  pm 
fraction:  diagonal  pattern  10-20%  of  sample,  solid  pattern 
^20%  of  sample.  Age  of  80 0 cm  sample  estimated  as  127,000 
years.  From  Sancetta  et  al.  (1973)  Figure  2. 

Figure  2.  Plots  of  f oraminif eral  index  Tw, S 018  in 
%o , and  foraminif eral  index  S versus  depth  for  Caribbean 
core  V12-122 . The  f oraniniferal  indices  are  defined  by  a 
set  of  transfer  functions  F3  and  scaled  to  be  unbiased 
estimates  of  winter  surface  water  temperature  (°c)  and 
average  surface  salinity  (%«),  respectively.  The  faunal 
index  curves  are  plots  of  a running  mean  of  three  samples. 
Letters  T through  Z designate  Ericson  zones  defined  on  the 
basis  of  presence  or  absence  of  £.  menardii . Letters  A through 
G designate  a correlation  with  the  climatic  cycles  of  Kukla 
(1970)  defined  from  the  soil  sequence  of  Central  Europe. 

From  Imbrie  et  al.  (1973)  Figure  3,  with  modifications. 

Figure  3.  Observed  values  of  winter  sea-surface  tem- 
perature versus  estimates  calculated  by  transfer  function 
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F13  for  191  North  Atlantic  sea-bed  samples  of  planktonic 
foraminifera . (a)  Scatter  diagram  with  80%  confidence  in- 

tervals indicated  by  dashed  lines.  (b)  Geographic  distri- 
bution of  the  residuals.  IA|  is  the  absolute  value  of  the 
difference  between  observed  and  estimated  temperatures. 

The  standard  error  is  1.165°C.  The  mean  deviation  is  indi- 
cated for  areas  with  more  than  two  points  in  which  more 
than  half  of  the  residuals  are  of  the  same  sign.  From  Kipp 
(1975),  Figure  28. 

Figure  4.  Sea-surface  temperatures  in  the  North 
Atlantic  Ocean.  (a)  Summer  temperatures  today  in  °C . (b) 

Summer  temperatures  18,000  years  ago  in  °C , as  estimated 
by  foraminiferal  and  coccolith  transfer  functions.  Solid 
contours  are  controlled  by  micropaleontological  data  in 
about  100  cores.  Sea-ict  margin  identified  from  paleonto- 
logic  and  sedimentologic  characteristics.  Note  that  this 
reconstruction  shows  substantial  cooling  in  northern  waters, 
with  little  or  no  change  in  the  subtropics.  From  McIntyre 


et  al.  (1975) . 
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